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1. Introduction 



The purpose of this note is to clarify the relation between the Hamiltonian and La- 
grangian approaches to nonlinear evolution equations. In particular we explain the least 
action principle and the Noether theorem in this context. The specific point of view, 
adapted from Souriau's book [8], has perhaps not been presented for nonlinear evolution 
equations. It seems that in the mathematics literature, as in [T],[l],[5], and in references 
given there, the Hamiltonian point of view is prevalent. In the physics literature, see for 
instance [3],[Z], the Lagrangian point of view rules with the symplectic structure largely 
neglected. In ^we present one possible mathematical reason for that. 



In this section we recall well known facts about the Hamiltonian structure of the nonlinear 
Schrodinger equation. The same point of view applies to other evolution equations, see for 
instance [1] and references given there. 

For simplicity we will consider the case of dimension one, and 



viewed as a real Hilbert space. The inner product and the symplectic form are given by 



Let H : V ^ W he a function, a Hamiltonian. The associated Hamiltonian vector field is a 
map Eh : V — ^ TV, which means that for a particular point u & V, we have (Eh)u ^ T^V. 
The vector field Eh is defined by the relation 



2. The Hamiltonian structure 



V = H\RX) C L2(M,C) , 



(2.1) 




(2.2) uj{v, {Eh)u) = duH{v) , 

where v G TuV , and duH : TuV — ^ M is defined by 

duH{v) = H{u + sv) . 



ds s=o 
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In the notation above 

(2.3) dH^iv) = {dH^,v) , {~H)u = -dH^. 

% 

If we take V = H^{M., C) with the symplectic form (12. ip . and 

H{u)= I -1^1^+1 



--dlu - \u\P-\ 



then we can compute 

duH{v) = Re y {{l/2)d^ud^v - \uf-^uv) 
= Rey {-{l/2)dlu- \u\''-^u)v. 

Thus, in view of (Q and (Q, 

{^h)u = \ (-^* 

The flow associated to this vector field (Hamiltonian flow) is 
(2.4) u=i^^^)^ = l.(^J_dlu-\ur\^ . 

3. The Lagrangian point of view and the Noether theorem 

According to [8] the following point of view towards dynamics goes back to Lagrange. 
We consider ^ 

y = -1/ X R = i/i(M,C) X R, 

and the following one form on V (we are rather informal here about dual spaces etc): 

(3.1) a^^,t){v,T) = K:{u,v)-H{u)T, {v,T) e T^^^t^V , {u,t)EV. 

Remark. The presence of the factor 1/2 in front of u in the definition of a is best 
understood using the finite dimensional analogy: if 2 = x + 2^, x, ^ G R"", then 

(3.2) a = - Im zdz — H{x,^)dt = -{^dx — xd^) — H{x,^)dt . 



We then define the differential of a: 



~ dcf , 

uj = da , 



that is 

25(«,t)((Wl,Ti), (W2,T2)) = Uj{Vi,V2) - dHu{vi)T2 + dHu{v2)T2 
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where we used the notation of ^ This calculation is easily understood using the analogy 
with 

d{^dx — xdS,)/2 = dC, A dx . 

Having uj makes V a presymplectic space in the sense that u has a kernel of dimension 
one. Here, the kernel is 

keTU(^,t) = {iv,T)eT(^^,t)V; ^ {v' ,T') e T(^^,t)V , T), K, T')) = o} . 

The following proposition replaces (12.41) with a condition related to cD = da: 
Proposition 1. The curve t ^ u{t) & V is a solution to 

(3.3) iut = -^dlu - \u\P~^u , 
if and only if 

(3.4) {u{t),l) ekeiUuit)- 
In other words, 

keruJu = ^C^u, 1) • 

Proof. We already know that (13. 3p is equivalent to (12. 4p . We then check that 
uj{{{Eh)u, 1), (v, T)) = uj{(Eh)u, v) - {dHu, T{Eh)u - v) 
= -{dHu, v) - {dHu, T(Eh)u - v) 
= T{dH^, {Eh)u) = T{dH^, (lA)rfi/„) = . 

□ 

A special case of Noether's Theorem (see [H (11.12)] for a more general version using 
the moment map) is now nicely given using this point of view: 

Proposition 2. Suppose that 

A{s) : {s,U)^U{s), sgM, UeV, 

is a one parameter group acting on V and preserving a: 

(3.5) A{sya = a, 

(here the pullback is given by f*ai^u,t){v,T) =^ afi^u,t){f*{'"iT))). Then 

(3.6) F{u,t) = a^u,t^(^^^A{s){u,t)\s=^ , {u,t)eV 
is conserved by the flow (13. 3p . 
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Proof. In the finite dimensional case we use Cartan's formula: if {d/ds)fs\s=Q = X (liere 
fs-V^ V), then at s = 0, 

= d{a{X)) + {da){X, •) . 

If we take = A{s) then the left hand side is and X = {d/ds)A{s){u,t)\s=o- The 
invariance of F is then equivalent to 

d{a{X)){u, 1) = 

but since 

d{a{X)) = -u{X,.), 

this follows from Proposition [1] The same argument applies formally in the case of evolution 
equations and can be easily verified. □ 

3.1. Standard group actions. The basic group action to consider are 

{u, t) I— s> {e~^^u, t) , (m, t) — s), t) , (m, t) {u,t — s) , 

and in each case we quickly see that A{s)*a = a. In the three cases we have 

{d/ds)A{s){u,t)\s=o = (-m,0) , {d/ds)A{s){u,t)\s=o = {-u^,0) , 
{d/ds)A{s){u,t)\s=o={0,-l), 

respectively, and the conserved quantities obtained using the formula (13.61) are easily seen 
to be 

J \u\'^dx , Im J u^udx , H{u) . 

A more interesting example is given by considering the Galilean invariance: 

A{s) (m, t) = {Ao{s, t)u, t) , Ao{s, t)u = e-''''l^+'*'u{» - st) . 
We first check that (13.51) holds. In fact, 

[A(s),](„,t)(t;,r) = {Ao{s,t)v + dt{Ao{s,t)u)T,T) = {Ao{s,t){v - + su,)T),T) , 

and hence 

(A(s)*a)(„,j)(t;,T) = aAis)(u,t){[Ms)*](u,t){^,T)) 

= uj{Ao{s, t)u, Ao{s, t){v - (is V2 + su^)T)) - 2H{Ao{s, t)u)T 
= a{u,t){v,T) , 

since 

Lj{Ao{s, t)u, Ao{s, t)v) = uj{u, v) , 

and 

H{Ao{s, t)u) = H{u) + s^(n, u) /4 + s(m, n^.) /2 

= H{u) + s'^uj{iu, u)/4 + suj{u, n^.)/2 . 
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We also see that 

-^A{s){u,t)\s=o = {ixu -tux,0) , 
as 

formula (13.61) gives 

Fiu, () = (!,„/ uM^ - / x| aNx = F(.. 0) = - / .1 al^d. = 

which of course corresponds to p = mq / 1 where 

p = Im I u^udx , q= — / x\u\'^dx , m = [ \u\'^dx , 



m _ 

are the momentum, position, and mass, respectively. 



3.2. Scaling. Let us now consider another group action preserving solutions of (13. 3p [p > 
1): 

(3.7) (s,M,t)i — > {Ao{s)u, s-H) , Ao{s)u{») = s^^u{s») . 
Then 

[Ais).]^u,t)i^,T)=Ais){v,T), 

and 

{A{sya)^u,t){^,T) = lu;{Ao{s)u,Ao{s)v) - H{Ao{s)u)s-^T = sl^{lu{u,v) - H{u)T) . 

That means that the form is preserved for p = 5. For p 7^ 5 we still preserve the kernel of 
uj = da which is consistent with (13. 7p preserving the solutions. 

To see the invariant quantity given by Noether's theorem (formula (13. 6p ) for p = 5 we 
compute 

^A{s){u,t)\s=i = {u/2 + xu^,-2t), 
ds 

and the conserved quantity is 

(3.8) F(m, t) = icu(M, u/2 + xu^) + 2tH{u) = -| Im j xu^u + 2tH{u) , 
which is a version of the virial identity, typically written 

(3.9) dt ( Im / xu^u] = AH{u) . 
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3.3. Case of p = 5. Here the scaling symmetry is part of a more general scaling property: 

see [B] for this and a recent study of the quintic NLS. 

Motivated by this, for g G 5-^2 (M) we define the standard action on M: 

, . at + b f a h \ ^tti,\ 



ct + d' \c d 

Then 

A{g) -.V^V, 

is given as follows 

A{g){u,t) = {A,{g)u,g~\t)) , Ao{g)u = (^?'(t))-ie-^^"W^'/(^^'W)u [{g'{t))'^x 

Since g'{t) = {ct + ig'{t))^ = {ct + d)~^ is well defined. 
The cases of 

' a \ fib 
l/a J ' \0 1 

correspond to scaling and translation with the invariant quantities already discussed. 
For 

/ X dcf / cos s — sin s 
9{'V = 

\ sm s cos s 

we obtain 

— {A{g{s)){u, t)\s=o = ((-t/2 + 2xV2)m - txu,, 1 + t') , 
ds 

so that the conserved quantity is 

F(u,t) = —- J x'^\u\'^dx + -tlra J xu^u — H{u)t^ — H{u) . 
Since H{u) is conserved and we also have (13. Sp we conclude that 
(3.10) 1 x^\u{x,t)\^dx= I x^\u{x,Q)\^dx + 2t\m I xu^u-AH{u)t^ 



which is again a version of the virial identity. This version of the virial identity is usually 
written 

Two time integrations, then substituting the identity 

dt I x^ImP (ix = 2 Im / xud^udx 
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evaluated at t = 0, give 

t + 4H(u)t'^ 





+ 2 Im / xudxU dx 




t=o J 



, which substituted 

i=0 



t=0 

Integrating (13. 9p from to t gives an expression for Im J xudxU dx 
here gives (13.101) . 

4. The least action principle 

To formulate the least action principle we need to define the Lagrangian. In the last 
section, although we took the Lagrangian point of view, we used the form a given by (13. ip . 
The Lagrangian, 

C ■ TV — ^ M, 



is defined as follows: 



/:(M,t,x,r) = «(„,j)(x,T), xeTuV. 



If 1 1-^ n is a curve in V we use a simplified notation 

(4.1) £(m) =^ a(„^f)(M, 1) . 
For the equation (13. 3p we obtain 

(4.2) C{u) = ]^Lu{u,u) - H{u) = -Um j ^tu-\j l«.f + ^ j l^r^ 

Action is more natural than considering Lagrangian. Let 7 be a curve in V . Then the 
action on 7 is defined as 

(4.3) S/= fa. 

When the curve is given by t t— ^ {u(t),t) we get, in the notation of (14. ip . 

Sj =^ J C{u)dt . 

The least action principle can be formulated as follows: 

Proposition 3. The curve 7 : s t— > (u(s),t(s)) is critical for if and only if 'y{s) G 
ker u^^s) ■ In other words 

(4.4) = <J==^ 7(s) G kera;^(^) . 
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Proof. We first give the proof in finite dimensions. Let 7,. be a smooth family of curves 
such that 7o = 0, and 7,. is equal to 7 outside of a compact subset, disjoint from d^y. Being 
stationary means that for any such family, 

d 



dr 



a 

7r 



= 0. 

r=0 



Let Fr be a smooth family of diffeomorphism such that, for r small, 7,. = Fr{'~f), and let 
X = {d/dr)Fr\r=o be a vector field defined on 7. Then, as in the proof of Proposition [H 
we use Cartan's formula: 
d 
dr 



a 

7r 



r=o dr 



— / -T^a 



r=0 



{da{X,m) + d{a{X))) 



[ u{X,.) + a{X)\9,= I u{X, 

J 'y J 'y 



'7 -^7 

since by the assumptions on 7^, X = near ^7. This means that 
d 



a 



= =^ 25^(,)(X^(,),7(s)) = VX V 

r=0 



dr 

which proves the proposition in finite dimensions. 



The same formal argument applies to evolution equation and in our case we check it by 
a sandard direct computation: 

S{u + 5u,t + 6t) = j {uj{u + Su,u + Su)/2 - H{u + 5u){i + Si))ds . 

Integrating by parts and neglecting higher order terms we obtain the first variation of S: 

SS = J {u{Su, u) — iduH(Su) + duH{u)5t)ds 

uo(u,t){{5u, St), {u,i))ds, 
and this vanishes for all 5u and 6t if and only if {u,t) G keTU(^u,t)- D 

5. Effective dynamics 

Suppose that M C is a submanifold which is presymplectic in the sense that 
(5.1) dimkeru5|^ = k , where k is constant on M. 

Then keTuj\j^ defines a foliation of M with leaves of dimension k. We note that the fact 
that du = and the formula for dp{X, Y, Z), 

Xp{Y, Z) - Fp(X, Z) + Zp{X, Y) - p([F, Z], X) + p([X, Z^Y) - p([X, F], Z) , 

show that the kercD satisfies the Frobenius integrability condition. 
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The method of collective coordinates for motion close to M is based on the following 
principle: 

Suppose that 7 is critical for S (for instance t u{t) which satisfies (13.41) 
or, equivalently, (12 ■4p ). Suppose also that 7 is close to M. Then it is close 
to a fixed leaf of the above foliation. 

Here is a trivial example to illustrate this. Let V = T*]R and H{x,^) = ^'^/2. Then 
suppose that M = = 0}. In that case dimkeTuj\j^ is 3. If 

7(t) = {{x + te,e),t), 

then it is close M, which the only leaf of the foliation. 

What one normally wants is (see [7] for examples from the physics literature and [2] for 
an implicit application of this principle in the mathematics literature): 

Let M satisfy (15.11) with k = 1. Suppose that 7 is critical for S. Suppose 
also that 7 is close to M. Then 7 is close to a 7^ C M which is critical for 

i5^Q, 7o C M. In other words, we restrict the Lagrangian to the submanifold 
and compute the action there. 

The simplest case is given by M = M x M with M symplectic, that is M for which u\m 
is nondegenerate. In that case the foliation is given by 

s {exp{sEH\J,s) , 

where Hp is the Hamilton vector field of a Hamiltonian p. This is very clear in finite 
dimensions since then, locally, 

M = {x,o--x"=e = o, x = ix',x"), e = (e',n, 

and 

^\m= (dC + H^'{x', 0, 0)dt) A {dx' - H^^x', 0, 0)dt) . 

However we may have situations in which u\m is degenerate yet keru;|MxR keeps fixed 
rank 1. That means that the Hamiltonian formalism is not applicable but the Lagrangian 
one is. Here is a simple example: 

M = {(xi, 0, ^2', ^2)}CV = T*R^ , H{x, 0=^1, 
uj\m = 26c?6 A dxi , dimker co'|A,/n{a=o} = 2 = dimM , 
= (2^2'^^2 + dt) A dxi dimkeru;|j^ = 1 . 

A complicated example from the physics literature comes from [21 §5]. If one considers 
M G V = H^(R,C) given by all sums of time- modulated solitons and time modulated 



10 



J. HOLMER AND M. ZWORSKI 




Figure 1. The plot of det{u\Mn{v=z=o}) for 7 = 0.1 (in the notation of [3]), 
(j) = 0, tp = 7r/4, < ?7 < 1.2 and 0.11 < a < 1.2. The two hnes afong which 
uj\Mn{v=z=o} is degenerate are clearly visible. The restriction to Z = V = 
is not essential since these two variables are essentially conjugate. 



nonlinear ground states (see [3], (5.1)]), 

u{x) = us{x; 7], Z, V, 0) + Unix; a, 0, ip) , 
(5.2) us{x; 7], Z, V, 0) = r] sech(r/x - Z)e^^^-^'^ , 

ud{x; a, 0, -0) =^ a sech (^ax + tanh~^ (~)) > 

then in the reduced six dimensional space described by {rj, Z, V, 0, a, ip) is not symplectic 
with respect to uj given by (12.11) . This can be checked by computing the determinant of a 
matrix corresponding to uj\m ^ see FiglH Despite that the method of collective coordinates 
is used by the authors in constructing an effective Lagrangian [31 (5.4)] and it is then used 
to obtained approximate equations of motion [31 (5.9)]. It is numerically shown to give a 
good agreement with the solution of the equation. 
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Finally we comment on the effective dynamics of solitons interacting with slowly varying 
potentials. In [5] we followed [1] and used a symplectic approach improving the results of [1] 
and [2] (same method apply to in that setting) by obtaining equations of motion without 
errors and obtaining a better accuracy of approximation by a moving soliton [h — > h"^, 
where h is slowness parameter of the potential). When attempting to reproduce these 
results using the Lagrangian formalism we could obtain the same equations of motions (we 
later learned that they were implicit in [7]), but could not obtain the h ^ improvement. 
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